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Abstract 

To estimate realistic cross sections in ultra peripheral heavy ion collisions we must remove 
effects of strong absorption. One method to eliminate these effects make use of a Glauber model 
calculation, where the nucleon-nucleon energy dependent cross sections at small impact parameter 
are suppressed. In another method we impose a geometrical cut on the minimal impact parameter 
of the nuclear collision (6 m , n > R\ + R2, where Ri is the radius of ion "i"). In this last case the 
effect of a possible nuclear radius dependence with the energy has not been considered in detail 
up to now. Here we introduce this effect showing that for final states with small invariant mass 
the effect is negligible. However when the final state has a relatively large invariant mass, e.g. an 
intermediate mass Higgs boson, the cross section can decrease up to 50%. 
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Collisions at relativistic heavy ion colliders like the Relativistic Heavy Ion Collider 
RHIC/Brookhaven and the Large Hadron Collider LHC/CERN (operating in its heavy ion 
mode) are mainly devoted to the search of a quark-gluon plasma in central nuclear reactions. 
In addition to this important feature of heavy-ion colliders, ultra peripheral collisions may 
give rise to a huge luminosity of photons opening the possibilities of studying two-photon 
and other interactions as reviewed in Refs.fl, 2, 3]. In the early papers on peripheral heavy 
ion collisions the effect of strong absorption was not taken into account. The separation 
of the strong interactions effects was solved by using impact parameter space methods in 
Refs. [4, 5, 6]. In order to obtain a truly peripheral photon-photon interaction one has 
to remove completely the central collisions, i.e we must enforce that in the cross section 
calculation the minimum impact parameter, b m i n , should be larger than Ri + R 2 , where 
Ri is the nuclear radius of the ion "i" [4]. The photon distributions can be described us- 
ing the equivalent-photon approximation (EPA) with the requirement of minimum impact 
parameter (or geometric cut) discussed above [3, 6]. 

The above method is not the only manner to avoid events where hadronic particle pro- 
duction overshadows the 7 — 7 interaction, i.e., events where the nuclei physically collide. 
An alternative is to use the Glauber model for heavy ion collisions [7]. It is a semiclassical 
model picturing the nuclei moving in a straight path along the collision direction, and gives 
the nucleus-nucleus interaction in terms of the interaction between the constituent nucleons 
and nuclear density distributions. If we write the cross section for the collision of two nucleus 
A and B as a function of the impact parameter (6) , the elastic (eZ) peripheral cross section 
will be given by 

a el = J d 2 b[l - exp (-ABa T AB (b)/2)] 2 , (1) 
where A and B are the nucleon numbers, o"o is the total nucleon-nucleon cross section and 

Tab® = J 70y^(g 2 )F B (g 2 )e^, (2) 

where Fa(b) are nuclear form factors. Eq.(l) and the form (2) for Tab(P) are valid only if 
one can neglect the finite range of the nuclear interaction. If at higher energies the total 
cross section increases both due to strength and due to the range the equation for Tab(P) 
should take this into account. The exponential factor in Eq.(l) is the one responsible for 
the suppression of the inelastic collisions. The o"o total nucleon-nucleon cross section that 
appears in Eq.(l) is known to be dependent on the energy. Actually the increase of hadron- 
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hadron total cross sections have been theoretically predicted many years ago [8] and these 
predictions have been accurately verified by experiment [9]. For instance, the proton-proton 
total cross section roughly double as we go from a few GeV up to the Tevatron energies. 

In ultra peripheral heavy ion collisions it is clear how this energy dependence of the 
cross section enters in the Glauber approximation. However the same is not true when 
we compute the cross sections with the EPA and the requirement of a minimum impact 
parameter. It seems that cross sections in very peripheral heavy ion collisions calculated 
within the Glauber method turn out to be slightly different from the ones computed with 
the geometric cut [10]. 

The nuclear radius certainly expands with the increase of the energy in the same way 
as the proton expands, and this expansion should be implemented in the geometrical cut 
calculation of peripheral heavy ion collisions. As far as we know this effect has not been 
discussed in detail in the literature, and it is the purpose of this work to introduce the energy 
dependence of the nuclear radius in the calculations of peripheral heavy ion collisions when 
the geometric cut method is used. 

In order to introduce the energy dependence of the nuclear radius in the calculations of 
peripheral heavy ion collisions we start discussing a standard computation of the photon 
distribution in the ion with the geometric cut method. The photon distribution in the 
nucleus can be described using the Weizsacker- Williams approximation (or EPA) in the 
impact parameter space. Denoting by F(x)dx the number of photons carrying a fraction 
between x and x + dx of the total momentum of a nucleus of charge Ze, we can define the 
two-photon luminosity through 



where r = s/s, s is the square of the center of mass (c.m.s.) system energy of the two 
photons and s of the ion-ion system. The total cross section of the process ZZ — > ZZX is 



where a(s) is the cross-section of the subprocess 77 — > X. There remains only to determine 
F(x). In the literature there are several approaches for doing so, and we choose the con- 
servative and more realistic photon distribution of Ref.[6]. Cahn and Jackson [6], using a 
prescription proposed by Baur [4], obtained a photon distribution which is not factorizable. 




(3) 




(4) 
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However, they were able to give a fit for the differential luminosity which is quite useful in 
practical calculations: 



where z = 2MRy/r, M is the nucleus mass, R its radius and £ (z) is given by 

^)=f>e- fe - (6) 
i=i 

which is a fit resulting from the numerical integration of the photon distribution, accurate to 
2% or better for 0.05 < z < 5.0, and where A l = 1.909, A 2 = 12.35, A 3 = 46.28, h = 2.566, 
b 2 = 4.948, and b 3 = 15.21. For z < 0.05 we use the expression (see Ref. [6]) 



Eq.(3) is written in a factorised form, which of course is valid only if one neglects the 
exclusion of central collisions into account. Therefore Eq.(3) is not the most general form, 
[3, 6], and the same is true for Eq.(5). The calculation assumes that the same radius R is used 
for both ions b min = 2R but also to have a cutoff for the individual impact parameter bi and 
b 2 (which either is necessary to eliminate final state interaction, or which takes into account 
the form factor effects, that is, the decrease of the charge inside the nucleus). Especially when 
looking, for instance, an intermediate mass Higgs boson production or other non-strongly 
interacting particles there is no reason to assume that the size of the individual cutoff radii 
for b\ and b 2 scales in the same way as b m i n . Therefore the calculation overestimates the 
dependence on R a bit. 

The condition for realistic peripheral collisions (b m i n > R\ + R 2 ) is present in the pho- 
ton distributions showed above. To obtain the above equations an elastic Gaussian form 
factor and an energy independent nuclear radius giving by R ion = 1.2 A 1 ^ fm have been 
used. A more accurate Woods-Saxon distribution for symmetric nuclei would produce some 
small deviations, but for our purposes the expressions for the luminosity described above 
are enough. However the expression for the nuclear radius is exactly the one we believe 
that should be changed by its energy dependent expression, and the problem is to have a 
phenomenologically sensible expression for the nuclear radius increase with the energy. 

In the heavy ions colliders nucleus like Au and Pb will collide with a great amount of 
energy, going from 200 GeV/nucleon (Au at RHIC) up to 5.5 TeV/nucleon (Pb at LHC), 




(5) 




(7) 
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and the ultra peripheral collisions can be computed with the help of the photon distribution 
described above. If the ion radius increase with the energy, the value corresponding to b m % n 
will also become greater, and consequently the cross section must decrease. This is easily 
seen in the many examples calculated in the literature where the cross section for a given 
process is concentrated at some moderate impact parameter and decreases when b increases. 
Of course the Lorentz factor is also important to determine this behavior. Therefore, if we 
introduce the energy dependence in the nuclear radius we expect lower rates for a given 
process than those obtained in the usual calculations, and this effect, even if it is small, 
could be important if we have a high precision measurement. 

The authors of Ref. [11] modelled the particle production process in ultrarelativistic 
heavy ion collisions in terms of an effective scalar field produced by the colliding objects, 
in their work they showed that the nuclear cross sections increase with the energy due to 
a logarithmic increase of the nuclear radius with the energy. We shall use this reference to 
obtain a relation between the nuclear radius and the incident energy that is the following: 

RUs)=1+2 ^ E+{S+l) ^ n (^y (8) 

R' = R P + R T ~ 2.4A 1 / 3 fm (R P (R T ) means projectile (target) and we assume Rp = 
Rt = Rion), is the energy of the projectile nucleus in the laboratory rest frame. The 
nuclear density for a nucleus A at distance x from its center is modelled by a Woods-Saxon 
distribution for symmetric nuclei, 



pws{x) = T ?Lr rrv ( 9 ) 

(1 + exp ^ j ' I) 

where d = 0.549 fm, and po can be found when the Wood-Saxon density is normalized by 
the condition / d 3 xp(x) = A. And e is equal to 



e = M z d 



16 R' (RtRpY 5 ' 7 ^ 
* 2 9 2 P% d? r 2 (*±i) 



(10) 



M z is the nuclear mass. The coupling constant g and the parameter for the mass spectrum 
5 were estimated in Ref. [11] and they are equal to 5 = —0.56 and g = 3.62 fm( 7+5 )/ 2 . 

The radius R 2 H (s) appearing in Eq.(8) at small energies gives a nuclear radius larger than 
Ri on , for this reason we have assumed the following normalization 

R 2 (s) RH{ S )R 2 ion , U N 



Ion 




R-ion 


R(a) 


Au 


0.2 


6.98 


7.29 


Ca 


7.2 


4.10 


4.75 


Pb 


5.5 


7.11 


7.61 



TABLE I: Values for R ion and R(s), Eq.(ll), in fm. The energies (y/s) are in TeV/nucleon. 

where Ri on = 1-2 A 1 / 3 fm. With this normalization factor we assure that when the ion energy 
is equal to its mass, the nuclear radius will be equal to Rion- It is the radius given by Eq.(ll) 
that should be used in Eqs.(5 - 7). Typical values for R ion and R(s) are showed in Table I. 

To show the effects of the nuclear radius dependence on the energy we computed produc- 
tion of leptons pairs (muons and taus) and resonances formed by the photon-photon fusion. 
In the resonance case we considered the r] c meson and an intermediate mass Higgs boson 
with a mass equal to 115 GeV. We computed the cross sections for two cases: In one the 
nuclear radius is energy independent (and equal to Ri on )- In the second case the radius 
obeys Eq.(ll). 

For an invariant mass of the photon pair above the threshold y/§ > 2m/, a lepton pair can 
be produced in two-photon collisions (77 — > l + l~) and the lowest order QED cross section 
for this subprocess is given by [3] 



^(77 - / + n = —a 



(3 -A 4 ) fl + P, 



Wi Vi -A 



(12) 



where A = \J 1 — Amf/s is the velocity of the pair in the 77 rest frame, m/ is the lepton 
mass, a/§ is the cm. system energy of the two photons and a is the fine-structure constant. 
Using this elementary cross section in Eq.(4) we obtained the rates shown in Table II. The 
calculation was performed for three different ions with different beam energies, the one of 
RHIC (Au) and the ones expected at LHC (Ca and Pb). The cross sections were integrated 
in a bin of energy equal to 1 < < 10 GeV. The third column of Table II shows the 
cross section computed with a constant nuclear radius and the fourth column the one with 
the energy dependent radius. For the three different ions the cross sections decrease when 
we consider the energy dependent radius described by the Eq.(ll). In all the cases the 
decrease is smaller than 10% and is negligible considering the theoretical and experimental 
uncertainties involved in the problem. 
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Ion 




(Jr. 

1 Hon 


a R(s) 


Ratio 


Au 


0.2 


2.127 


1.947 


1.09 


Ca 


7.2 


0.643 


0.588 


1.09 


Pb 


5.5 


106.4 


101.3 


1.05 



TABLE II: Cross sections of the process ZZ77 — ► ZZ[i + fi~. The cross sections u Ri (ctr( s )) given 
in the third (fourth) column are the ones computed with the energy independent (dependent) 
radius. The last column shows the ratio between the third and fourth columns. The cross sections 
are in mbarn and the energies (y/s) are in TeV/nucleon. 



Ion y/s a Riori o R(s) Ratio 

Au 0.2 6.972 x 10~ 4 5.727 x 10~ 4 1.22 

Ca 7.2 5.176 x 10~ 3 4.604 x 10" 3 1.12 

Pb 5.5 0.759 0.718 1.05 

TABLE III: The same as in Table II, but for the subprocess 77 — ► t + t~. 

In Table III it is possible to see the results when the subprocess analyzed is 77 — > t + t~ 
with 2m T < < lOGeV. The general behavior of the r pair production cross sections is 
very similar to the one observed previously in Table II. Of course, the cross sections for 
producing tau pairs are smaller. However the collision of Au-Au and Ca-Ca are now more 
sensitive to the energy dependence of the nuclear radius, producing an effect larger than 
10%. The rates for tau pairs production in Pb collision with a cm. energy equal to 5.5 
TeV/nucleon, with and without energy dependence in the ion radius are not so different. 
As we shall discuss later the larger cut that we perform in the impact parameter when we 
consider the energy dependent radius removes photons of larger energy. Therefore for final 
states with larger invariant masses we may expect a larger effect. 

Let us now consider the case of heavy resonances. To estimate the production of one 
resonance R formed by a photon-photon fusion in peripheral heavy ion collisions we use the 
following elementary cross section in Eq.(4), 

(7(77 - R) = ^-/(R - 77)5 (V - ^pj (13) 
where M R is the resonance mass and F(R — > 77) its decay width into two photons. In 
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Ion y/H a Rxon a R{s) Ratio 

Au 0.2 2.147 x 10~ 3 1.846 x 10" 3 1.16 

Ca 7.2 2.897 x 10~ 3 2.614 x 10~ 3 1.11 

Pb 5.5 0.437 0.413 1.06 

TABLE IV: The same as in Table II, but for the subprocess 77 — > rj c . 

Ion yfi a Rion a R{s) Ratio 

Ca 7.2 9.970 x lO" 10 6.789 x lO" 10 1.47 
Pb 5.5 1.854 x 10~ 8 1.387 x 10~ 8 1.34 

TABLE V: The same as in Table II, but for the subprocess 77 — > H. 

Table IV we show the results obtained for two-photon production of i] c in peripheral heavy 
ion collisions with M Vc = 2.979 GeV and T(r] c — > 77) = 6.6 keV. The ratio of the cross 
sections considering the two scenarios are 1.16 and 1.11 for Au and Ca ions, respectively, 
and 1.06 for the Pb ion. Finally, in Table V it can be observed the values corresponding to 
the subprocess 77 — > H with M H = 115 GeV, where we used the Higgs boson two-photon 
decay width found in Ref. [12]. We do not show the result for RHIC energies because it is 
too small. The values of Table V indicate that the production cross sections for both ions 
are strongly affected by the inclusion of a radius described by Eq.(ll). In the case of Ca 
collision with a cm. energy of 7.2 TeV/nucleon the cross sections decrease nearly to half of 
the value obtained in the case of a energy independent radius. The situation is less drastic 
for the Pb ion with y/s = 5.5 TeV/nucleon, but the ratio is still large (= 1.34). This is the 
only situation that we investigated where the Pb collision is clearly sensitive to the use of 
Eq.(ll) (or to the energy dependence of the nuclear radius). 

The fact that a sharp cutoff in impact parameter space at b m i n should be replaced by a 
smooth one was already discussed in Ref. [2] . Comparing the Glauber model calculation with 
the one with a sharp cutoff we could expect significant deviations present at the upper end 
of the invariant mass distribution. Looking at the photon luminosity we see that only the 
smallest impact parameter contribute significantly to the events with large invariant masses. 
Imposing the cut on b min but now with the radius described by Eq.(ll) we obtain a more 
realistic calculation of the very peripheral heavy ion collisions. 
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In conclusion, we discussed the two different ways to compute cross sections for ultra 
peripheral heavy ion collisions. In the Glauber method it is quite clear how the increase with 
the energy of the nucleon-nucleon cross section enters in the calculation. In the calculation 
with the geometrical cut imposed on the impact parameter, the nucleon, as well as the 
nuclear, radius expansion with the energy was not introduced up to now. It was noticed 
in Ref. [10] that there was a difference between the two methods. The difference was small 
and had some dependence on the invariant mass of the final states. The work of Ref. [1 1] 
prescribe a very precise way to introduce the nuclear radius dependence with the energy. 

We believe that the estimative of the cross sections in ultra peripheral collisions with 
the geometrical cut method just changing the radius independent of the energy by the one 
dependent of the energy will give realistic predictions for any invariant mass of the final 
state. The effect is of order of 50% for an intermediate mass Higgs boson. Turning the 
problem the other way around we may also say that if the ultra peripheral collisions are 
measured with high precision, we may have a new way to study the increase of the nuclear 
radius with the energy. To do so we just have to measure the cross sections for very known 
final states with small and large invariant masses with high precision, there should be a 
decrease of the cross sections as a function of the invariant mass as we go to larger and 
larger energies. 

Note added: Some comments on the effects discussed in this work were also made by 
S. Klein and J. Nystrand in [13], where the Fig. 3 gives the reduction in gamma-gamma 
luminosity (for gold at RHIC) for a Glauber calculation of hadronic interactions compared 
to the one with geometrical cut. 
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